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We formulate a field theory capable of describing the en- 
tanglement in a canonical ensemble of polymers in terms 
of Feynman diagrams. 



I. INTRODUCTION 

The discovery of the intimate relationship between the 
statistical mechanics of long polymer molecules and cer- 
tain field theories describing critical systems has been 
crucial in achieving the present understanding of the 
physics of unentangled polymer chains [|l| . 

Under certain experimental conditions, polymers form 
knotted configurations which remain stable in time. 
In this case, the topological relationships among the 
molecules become relevant. The principal theoretical 
tools for investigationg these were set up by S. Edwards 
1^ and Brereton and Shaw ||^, and formulated by one of 
the authors (HK) with the help of the topological field 
theory of Chern and Simons Q . While the tools are avail- 
able, a satisfactory theoretical treatment of an ensemble 
of topologically linked polymers is still missing, in spite 
of much work [|5|~p^ . 

Definite progress has recently been achieved by the 
present authors by mapping the statistical me- 

chanics of two fluctuating polymers to a topological 
Ginzburg-Landau model, and applying Feynman dia- 
gram techniques to calculate the average of the square 
winding number. In the present work we extend the field 
theoretic formulation of [|ll|,0 to a system of N poly- 
mers. 

Excluded volume forces, although physically impor- 
tant, have so far been neglected, and will be ignored also 
in this work. 

We start in Section || from the path integral descrip- 
tion of polymers in which these are viewed as Brownian 
trajectories. The topological constraints are imposed us- 
ing the simplest link invariant; the Gauss link integral 
for pairs of trajectories. The topological interactions are 
complicated, and several nontrivial steps are necessary to 



arrive at a useful field theoretic description of the system. 
For this purpose we introduce in Section III a set of auxil- 
iary abelian Chern-Simons fields. These allow us to make 
the theory local and to convert the polymer path integral 
to a Markoffian form. Actually, there are several possible 
abelian Chern-Simons field theories which could accom- 
plish this task, differing from each other by the number 
of fields. However, we prove that these are equivalent 
after exploiting the field equations and the freedom of 
performing linear transformations of the fields. 

The arbitrariness in choosing the auxiliary topological 
field theory is used in Section VI to overcome a techni- 



*fferr ari@science.unitn.it 
tkleinert@physik.fu-berlin.de 
*lazi@tn. infn.it 



cal problem which was absent in the two polymer case: 
the coupling constants of the interactions between Chern- 
Simons fields and polymer trajectories are related to the 
parameters for the topological constraints by non-linear 
algebraic equations. The latter are too complicated to be 
solved analytically apart from particular cases, in which 
the parametrization simplifies considerably. A particular 
choice of the Chern-Simons fields removes this problem. 

With these methods we convert the polymer path in- 
tegral to a Markoffian form and it becomes possible to 
complete the mapping of the polymer problem into a field 
theoretical model of topological entanglement. This is 
achieved in Section^ by exploiting the method of replica 

0. 



II. PATH INTEGRAL APPROACH TO 
TOPOLOGICAL POLYMERS 

Let Pi, . . . ,Pn be a set of topologically linked poly- 
mers of lengths Li, . . . Ln respectively. In order to keep 
our treatment as general as possible we consider both, 
open and closed chains of polymers. Strictly speaking, 
the entanglement of open chains is not really topological. 
Over a long time scale, it is always possible to disentangle 
open polymers, and the topology of the system is not con- 
served. However, if we restrict ourselves to intermediate 
time scales short compared to that of complete molecu- 
lar rearrangements, open chains will behave almost like 
closed chains. In this approximate sense we shall be able 
to apply our topological methods also to open chains. 
Moreover, the study of the difference between the be- 
haviors of open and closed chains will make it possible 
to gain some insight in the nature of topological inter- 
actions [pf. Thus, we shall start our considerations with 
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open polymers Pi with i = 1, . . . , whose end points He 
at x\ y'. For coinciding end points = y*, they become 
closed polymers running through x\ 

The relevant quantity to describe the statistical me- 
chanics of the polymers is the configurational probability 
G{„i}(x, y; L), which measures the probability to find the 
polymers Pi with end points at x' , y' in a given topolog- 
ical configuration {m}. To compactify the notations, we 
have collected the set of all end points in iV-dimensional 
multi- vectors x = (x^, . . . , x^), y = (y^, . . . y^), with an 
analogous vector notation for the lenghts of the polymers: 
L = {Li,...,Ln). To distinguish different topological 
configurations, we shall use the Gauss link invariant: 
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ds, 



dsjx(si) 



, x'(g»)-xJ(g,) 



defined for each pair Pi, Pj of polymers with i ^ j. 

The integral (|l]) is well-defined also for open trajecto- 
ries, but it becomes a real topological invariant only for 
a pair of closed polymers with x' ^ y*, x' ^ y'. In the 
latter case it counts how many times Pi winds up around 

Following the approach of Edwards , the configura- 
tional probability can be expressed as a path integral over 
all possible configurations x*(si) with < Si < Li and pe- 
riodic boundary conditions x'(0) = y* and x'(Lj) = x': 



G 



{m} 



Px'(si) • ■ • / Px^(sAr)e-(-^°+-^=-) 

iV-l TV 

x]lI[^(x{P^,P,)-m^,), (2) 

i=l i=2 

j>i 

where is the euclidean action of a random-chain 
. 3 ^ 



c' ^ds, 



(3) 







and 



(4) 



the steric repulsion between the chain elements. The 
parameter a denotes the length of the chain elements and 
v^^ represents an N xN matrix of coupling constants with 
the dimension of a volume. 

The (5-functions in the integrand of Eq. (||) enforce 
the topological constraints that the pairs of chains Pi , Pj 
wind around each other a number of times my . Since we 
are describing open polymers up to this point, these num- 
bers are continuous. Only for closed polymers will they 
become integer numbers. Then the Dirac (S-functions in 
Eq. (0) have to be replaced by Kronecker S's. 



In the following, it will be convenient to introduce an 
auxiliary probability G'{a}(x, y; L), from which the origi- 
nal G{m}(x, y; L) is obtained by a Fourier transformation 
with respect to the topological numbers: 

G{„}(x,y;L)=/ J] J] ^e-^-™-G{,}(x, y; i) . 



= 1 j = 2 
3> 



(5) 



The auxiliary probability has the advantage that its path 
integral representation 



G{A}(x,y;Lj = J ^ V^\s^) ■ ■ ■ j ^ Px^ 



(1) 



-(A)+-4„+Aop) 



(6) 



accounts for the topological constraints among the poly- 
mers by a source like term: 



N 



(7) 



i=l j=2 



Note that if in a formulation for closed polymers, where 
uriij are integer numbers and the Dirac (5-functions in 
Eq. ^ are Kronecker symbols, the Fourier variable Ay- 
would be cyclic with a range \ij G (0,27r). Then we 
would use angular variables Lpij rather than Xij , and with 
Eq. d) as 



G'{m}(x; L) 



nn 

i=l 3=2 



e-^'^-"-G{;,}(x;L). (8) 



where 



G 



^}(x;i) = lim G{„} 

x^y 



(9) 



Returning to (^) , we see that it has the form of a path in- 
tegral over the trajectories of a system of N particles per- 
forming a random walk. Thus we may exploit the duality 
between particles and fields, which is valid in statistical 
mechanics as well as quantum mechanics, and express 
the path integral in terms of fields. Although the gen- 
eral techniques are well-known |Q , the specific task here 
is complicated by the presence of the topological term 
(^), which is non-Markoffian and rather complicated. We 
shall solve this problem by introducing auxiliary fields, 
which allow us to rewrite the right-hand side of in a 
more tractable form. 



III. AUXILIARY FIELDS AND THE 
DECOUPLING OF TRAJECTORIES 

Before coming to this wc first reformulate the excluded 
volume interaction Av in Eq. (0) is a standard way, since 
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its non-Markoffian character. The trajectories are cou- 
pled to each other by the two-body potential in Eq. (jj). 
The interactions can be disentangled by introducing N 
real scalar fields (j)i . . .(pN with the euclidean action 



(10) 



It is obviously possible to rewrite the exponential of —Ay 
in ^ with the help of the following identity M: 



y • ■ •P0Are-^* J]^ X cxp -i 



dsiipi (x*(si)) 



(11) 



On the right-hand side, each trajectory x'(si) interacts 
only with a single field 0i(x), so that the polymers move 
in individual random fields. Thus, the contribution of 
the excluded volume forces is converted to a Markoffian 
form. 



IV. TOPOLOGICAL INTERACTIONS 

Let us now turn to the topological interactions, where 
we search for auxiliary fields to simplify .Atop in the path 
integral . Such auxiliary fields are provided by Chern- 
Simons theories Let A"^ with a = 1, . . . , TV' be a set 
of N' Chern-Simons fields with euclidean spatial indices 

= 1,2,3. They will allow us to write the identity 



JVAe 



-iAcs 



where is a the Chern-Simons action 



N' 



(fixA" • fV X A'^l 



a, 0=1 • 



and J* are currents 



J^(x) 



ds,x*(s,)^'" (x-x^(s.O) 



(12) 



(13) 



(14) 



The measure of functional integration DA is short for 

the product Ha^i I'A". We assume gap to be a suitable 
N' X N' -symmetric matrix to be specified later, which 
possesses an inverse (g~^)"'', while hia is a. N x N' ma- 
trix. 

The right-hand side of Eq. ( p^ ) has the described 
Markoffian form. To calculate it explicitly, we quantize 
the Chern-Simons fields A'^ in the Lorentz gauge, where 
they are completely transverse, and have the correlation 
functions 



G"5(x,y)^(A"(x)A^(y)) 



l\al3 



{x - vY 

After some calculations we find 



(15) 



/PAe- 



■iAcs 



exp 



N 



(16) 



The right-hand side is equal to that of Eq. ( p_2[ ) with the 
topological action (|^) if wc satisfy the equation 



(17) 



Since the matrix Ay has vanishing diagonal elements, 
this condition automatically eliminates the appearance 
of non-topological terms of the form x{PiTPi) with i = 

At this point, it is convenient to introduce the linear 
combinations of the Chern-Simons fields 



N' 



such that Eq. ( [T^ ) becomes 



JVAe- 



-iAcs 



(18) 



(19) 



The absence of the non-topological contributions 
x{Pi,Pi) in (^6|) is ensured by requiring purely off- 
diagonal correlation with 



(C;(x)C:(y))=0, z = l,...,7V. 



(20) 



The linking of a polymer Pi with a different polymer Pj 
is described by the off-diagonal correlation functions, the 
related propagator should be different from zero: 



(C;(x)C7^(y)) ^0 i^j 



(21) 



This condition allows us to conclude that in the most gen- 
eral situation in which all elements Xij are non- vanishing, 
the number N' of Chern-Simons fields should be at least 
equal to N . 

Indeed let us assume the contrary: N' < N . In this 
case we can always rearrange the indices such that the 
first N' fields C'^(x) with cr = 1, . . . , A^' are independent. 
The remaining fields 



C^'+^(x) 



r : 



1, 



(22) 



should then be linear combinations of the first N' fields. 



N' 



(7;(x)=.^s;c;(x) 



(23) 



CT = 1 
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From the property (kG) , we find the correlation functions 



N 



N' 



(C;(x)C:(y)) = ^Wa'K{^)C:\y))=Q (24) 

a,a' — 1 

By hypothesis, however, we have 

(25) 



(C;(x)C,-'(y)) cx A,,-e^.p|^-^ ^ 0. 



As a consequence, the most general solution of Eq. (p4) 



if (TT^CT-, 



(26) 



where ct is a fixed integer with 1 < a < N' . Therefore, 
the fields C'^(x) and C°'(x) coincide apart from an irrel- 
evant factor Sj. In this way, if A^' < A^, we obtain for 
N' < N the following contradiction: 

= {c;{^)c:{y)) - 4{c;{^)c:{y)) 

= (C';(x)Cf(y)} = (Cf+^(x)C;(y)) 



{x-y)P 



(27) 



In the opposite case of A^' > A, it is always possible to 
reduce the number of Chern-Simons fields to A^. To show 
this we consider the action 



N 



(28) 



appearing in Eq. (p^. In the above action we perform a 
change of variables in which a number A^ of fields A" is 
expressed as a linear combination of the fields C"s and of 
the remaining A"'s. Without loss of generality, we may 
suppose ( [l^ ) to be invertible, with the solutions 



N / N' 

i=i 



(29) 



a=N+l 



Substituting this into (E8|), we find 



N 



N' N N' 
a=N+l 1=1 a.,f3=N+l 



Ala 



4" 



N' 

E 



(31) 



It is easy to see that the redundant fields A'" with a = 
A + 1, . . . , A' can now be integrated out is Eq. (p^, so 
that we arrive at a Chern-Simons field theory with A 
fields Ci,...,C^. 

This result could be expected from the fact that only 
the fields have sources in the action (p8|), while the 
remaining fields are free, and may be eliminated via the 
equations of motion. As a consequence, the freedom 
in choosing the number of Chern-Simons fields is only 
apparent, because any Chern-Simons field theory with 
N' > N substituting Eq. ( |l9| ) is equivalent to a Chern- 
Simons field theory with A^ fields only. In principle, there 
is still some arbitrariness in the choice of the matrix el- 
ements g"^ and hia once A^' has been fixed, since their 
values are only constrained by ([T^). However, this ar- 
bitrariness reflects merely the possibility of performing 
linear transformations of the ^'^'s, and it is thus irrele- 
vant. 

To conclude this section, we compute the denominator 
in the right hand side of ([l9|). For this purpose we note 
that since gafs is a N' x N' symmetric matrix, it can 
always be expressed as follows: 



9af3 



N' 

E 

/3=1 



(32) 



where 77^/3 is again a, N' x N' symmetric matrix. Per- 
forming in Eq. (13) the substitution: 



N' 



the Chern-Simons action becomes 



(33) 



N' 



i.e., the dependence on ga/3 disappears. Thus: 



J DAe-^-^^^(^°) = [det(5a;3)] 



1/2, 



(34) 



(35) 



where c — J Vae '"^"^^ ^ is an irrelevant constant fac- 
tor. 



(30) 



where the constant coefficients M ^ , Nai and Oaf3 are 
functions of the matrix elements g"^ and hia ■ The mixed 
terms in ^^g, which are proportional to Nai, are elimi- 
nated by introducing the new field variables 



V. FIELD THEORY 

In the previous section we have seen that the path 
integral over the polymer trajectories can be converted 
to a Markoffian form via auxiliary fields. In rewriting the 
topological interactions there is some freedom in choosing 
the auxiliary fields by varying their number. Also the 
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parameters gap and hai are not completely fixed by the 
system of equations 

On the other hand, it has been shown that all abelian 
Chern-Simons field theories for which the relevant iden- 
tity ( p^ is satisfied are equivalent after exploiting the 
equations of motion and performing linear transforma- 
tions of the vector fields. 

The simplest way to ensure (|9|) is to choose 



where the currents J^(x) have been already defined in 
Eq. (|l|) and the fields C* of Eq. (|l|) have the following 
explicit expressions: 



N-l 



J2 ^N^-^' 



= K^AyA^' + B\ i = 2,...,N -1. 



(41) 



(36) 



for i,j — 1,...,N. This choice, however, makes the 
Fourier transformation (^) from the auxiliary probability 
G{Ai,(x, y;L) to the original configurational probability 
G{m}(x,y;L) too difficult for an analytic treatment. As 
a matter of fact, starting from Eq. (pq), we find with the 
help of (11: 



VA 



The factor c^^ in Eq. (|4^) is an irrelevant constant inde- 
pendent of \ij . 

Using Eq. (^Tj) and (^0|) in the expression of the aux- 
iliary probability (^ we obtain the path integral repre- 
sentation 



N 



p = [(detAy)i/2c 



N N 



X exp <^ j d^x ^ ^ Ay 



1=1 i = i 



An 



dx''(s.)A^(x*(sO) 



G{A} (x,y;L) {J[G{^\r;U\cl),,C')j, 
where 

G{x\r;U\cj),,C') = / Px'(s,) 

Jy' 

X e"[^o''" 'i''.-C*,(x-(^-.))+/o''''i^-.i*(^.)-C'(x'(s.))] 



(42) 



(37) and 



The denominator in the right hand side contains the term 

1 /2 

[det(Aij)] , which complicates the integration over the 
parameters Xij in (|^). 

Moreover, in the ansatz ( |36| ) the requirement that the 
matrix gij should be invertible cannot be guaranteed for 
all possible matrices Xij . 

The situation does not improve if we choose 



/:^.(X*(S,)) = —x"^{s,)+iUx\■S^)) 



(43) 



(44) 



N' 



9i] 



47r 



^2 
47r 



■mo- 



rn 



In Eq. (|42|), the expression in the expectation symbol 
( ) must be averaged with respect to the fields <j)i with 
i = 1, . . . , iV and the Chern-Simons fields A^ , with 
j = l,...,7V-l. 

The path integral in Eq. ( ^3| ) describes a Markofhan 
random walk of a particle immersed in an electromag- 
netic field (C% In analogy with the evolution kernel 
of a particle in quantum mechanics, G (x*, y% Li\(f)i^ 
satisfies a Schrodinger-like equation 



Here the elements 7]ij have a complicated dependence on 
the variables Ay , via the algebraic equations (|T7|). 

To solve these difficulties, we exploit the freedom of en- 
larging the number of topological vector fields. The sim- 
plest Chern-Simons field theory for our purpose contains 
N' = 2{N - 1) fields A\...A^^^ and Bi,...,B^-i. 
The action Acs is given by: 



dL, 6 ' 



where 



D,; = V + 



(45) 
(46) 



N-l 



i V — J. li 

i—l 

Equation ( |l9|) becomes now 



(39) 



(40) 



It is now convenient to consider its Laplace-transformed 
in the length parameter Li: 

G'{A}(x,y; fl) 

dLi---dLN e-Sr=i'^'^'G{A} (x,y;i) . (47) 

The parameters fii in Boltzmann-like factors control the 
growth of the polymers. Applying the Laplace transfor- 
mations to both sides of Eq. (E2|) , we find 
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N 



(48) 



where G (x\ y*; C*) is the Laplace transformed 
correlation function (^), obeying the stationary equa- 
tion 
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G(x\y';^,|0„C») =<5(x'^-y'). (49) 



Equation (49) can be solved in terms of fluctuating poly- 
mer fields V* , V'i as 0] 

G(x^y';^,|<^„C') ^^J Vi,,Vr^ 
xV'.(x')^*(y')e-^-'['^*''^^l (50) 
where ^poi is the polymer action 

A^o\['^*,i^^]^ j d?x ^|DVip + (/i, +i(/),)|V',f (5f) 

and Zi is the associated partition function 

j V^iV^*e--^'^'"^'^''''"^ (52) 

The integrations over the auxiliary fields 0; is compli- 
cated by the presence of the factor in Eq. (pO[), which 
makes them non-gaussian. This problem can be solved 
exploiting the method or replicas. 

VI. REPLICA FORMULATION 



For each pair of fields ipi, ^* we introduce a set of ni 
replica field ip^^ , V'j*"' with = f , . . . , n^. It is convenient 
to group the replica fields in n,jtuplets 

vi/. = (^^...>r), ** = (v*\...,V'r)- 

With these fields the correlation function (^0|) may be 
rewritten as follows: 



G(x\y^M,|0„e) 

= lim / P*,;P**Vr(x')^r"(y') e 



where we have set 

v^,v^* = / n ^c'^V'r"' 

ai = l 

and defined the replica field action 

■rii 



= -|D,V^,|2 + (^,+z0,)|5',|2 
D 



(53) 



(54) 



(55) 



The index in ( p3[ ) is a fixed replica index chosen arbi- 
trarily in the range I < < n^. 

The limit of zero replica number in ( p3| ) is performed 
by an analytic extrapolation. The path integral on the 
right hand side is calculated for integer values of rii and 
the result is then extrapolated analytically to the point 
Ui = 0. 

Combining everything, the auxiliary probability (|4S 
has the functional integral representation 



G{A}(x,y;M)= lim [ VAVB 

' N 



N 



(56) 



Note that in this expression the integration over the aux- 
iliary fields has become Gaussian. After a suitable 
rescaling of the fields -0°' , i we give the final result in 
the following form: 



G{A}(x,y;A^)= lim [ VAVB 



N 



i=l 



X vr(x)Vf-(y^)]e-^- (57) 



where 



AT 

Act = Ucs '^''^ [^^* (-°^ + 



and 



(58) 



(59) 



With respect to Eq. (^, the fields Vf", '0r' have been 
rescaled by a factor so that they acquire the 

canonical dimension = [ip*'^'] = 1/2 of usual scalar 
fields. 

Following [Q, we have introduced a mass parameter 
M and a parameter playing a similar role as the Planck 
constant in quantum mechanical path integrals h = ^.a/i- 
The value of M has been fixed with respect to the step 
length a by requiring the condition fi — 1. 

In this way, the action (|5^) becomes that of a stan- 
dard field theory with unit Planck constant. Moreover, 
Atot is a quadratic form with respect to the parameters 
\ij , so that the inverse Fourier transformations leading to 
the original configurational probability can be performed 
after a diagrammatic evaluation of the correlation func- 
tions (|5^). 

The actual calculations are left to a future publication. 
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